Running coupling effects for the singlet structure function g\ at small x 
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, The running of the QCD coupling is incorporated into the infrared evolution equations for the 

f*"^ ■ flavour structure function gi . The explicit expressions for g\ including the total resummation of the 

double-logarithmic contributions and accounting for the running coupling are obtained. We predict 
CNJ , that asymptotically g\ ~ x~ As , with the intercept As = 0.86, which is more than twice larger than 

the non-singlet intercept Ajvs = 0.4. The impact of the initial quark (5q) and gluon (8g) densities 
on the sign of gi at x <C 1 is discussed and explicit expressions relating 8q and 8g are obtained. 
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I. INTRODUCTION 



The flavour singlet structure function gi has been the object of intensive theoretical investigations. First g\ was 
calculated in Refs/0,0, where the LO DGLAP evolution equations 0,0 were used. Since that, the DGLAP approach 
has become the standard instrument for the theoretical description of the singlet and non-singlet components of g\ 
and provides good agreement between experimental and theoretical results (see e.g. Ref. and the recent review 
0). On the other hand, despite this agreement, it is known that from a theoretical point of view the DGLAP 
equations are not supposed to work well at x <C 1 because the expressions for the anomalous dimensions and the 
coefficient functions in these equations account only for a finite part of the NLO (~ l/w Tl+1 ), with n = 1, 2 (see jg). 
In particular, this means that all double-logarithmic (DL) contributions ~ (a s m 2 (l/;r)) Tl , n > 1 are neglected in the 
DGLAP expressions for g\. As a matter of fact, these contributions become very important in the small x region 
and should be accounted for to all orders in a s . Such total resummation was done in Refs. for the non-singlet 
component, g^ s , of gi and in Ref. for the singlet one. These calculations were done in the double-logarithmic 
approximation (DLA) and it was shown that g\ has the power-like (Regge) behaviour ~ x~ A when x — > 0. However, 
k> the QCD coupling in Refs. 00 is kept fixed, whereas, as well known, the running coupling effects are relevant. As 
the DL intercepts A of g\ obtained in are proportional to y/a^ fixed at an unknown scale, it makes the results 
of Refs. 00 unclear and not suitable for practical use. A p arametrisation for fixed a s — a s (Q 2 ) for both the singlet 
and non-singlet component was sugg ested in Refs Q3- [lj. In the DGLAP framework, in a general ladder rung a s 
depends on the transverse momentum k± of the ladder parton (a quark or a gluon) , with /i 2 <k\ < Q 2 , where /i 2 is 
the starting point of the Q' 2 -evolution. The arguments in favour of such a dependence were given in Ref. 9J. However, 
in Ref. |l5fit was shown that the arguments of Ref. are valid only when x is large (x ~ 1), while on the contrary 
in the small- a: region, in a general ladder rung i a s depends on the virtuality of the "horizontal" gluon, (k{ — /c^-i) 2 . 
This dependence was used in Refs. 01 for studying g^ s at small x, with running a s accounted for. Our prediction 
for the non-singlet intercept was confirmed by phenomenological analyses of the experimental data in Refs. [Tsj . 

In the present paper we apply the same arguments of Ref. |15| in order to account for the running a s effects for 
the description of the small- a; behaviour of the flavour singlet component of g±. The paper is organised as follows: 
In Sect. 2 we construct and solve the system of the infrared evolution eqs. for g\. In Sect. 3 we obtain the explicit 
expressions for the anomalous dimensions of g\. These expressions account for all DL contributions and also for 
running a s . The intercept of g\ is calculated in Sect. 4. In Sect. 5, we specify the general solutions obtained in Sect. 2 
through the phenomenological quark and gluon inputs and estimate the sign of g\ at x -C 1. Then we calculate the 
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perturbative contributions to these inputs. Finally, Sect. 6 contains our concluding remarks. 



II. IREE FOR THE STRUCTURE FUNCTION gi 

In order to obtain g\ with all DL contributions accounted for at small x, we cannot use the DGLAP equations. 
Instead of that, we construct for it a set of infrared evolution equations (IREE), i.e. equations for the evolution with 
respect to the infrared cut-off /i in the transverse momentum space. The cut-off fi is introduced as the starting point of 
the evolution with respect to both Q 2 and x. In contrast to the DGLAP equations where only Q 2 -evolution is studied, 
the IREE are two-dimensional. In order to derive these IREE, it is convenient to operate with the spin-dependent 
invariant amplitude M q which can be extracted from the forward Compton scattering amplitude by using the 
projection operator it^vXplxVplvi- Through this paper we use the standard notations, so that q\ is the momentum 
of the off-shell photon (— q 2 = Q 2 > p 2 ) and p^ is the momentum of the (nearly) on-shcll quark. When M q (s, Q 2 ) is 
obtained, the polarized quark contribution to g%, g q , can be easily found: 



g q (x,Q 2 ) = -%M q (s,Q 2 ) (1) 

where s is the standard Mandelstam variable, s ss 2pq. The subscript q in the rhs of Eq. JIJ means that the off-shell 
photon is scattered by the quark. It turns out that the IREE for M q (x,Q 2 ) involve the spin-dependent invariant 
amplitude, M g (s,Q 2 ) of the forward Compton scattering where the off-shell photon is scattered by a nearly on-shell 
gluon with momentum p. Similarly to Eq. the polarized gluon contribution, g g is related to M g as follows: 

g g (x,Q 2 ) = -%M g (s,Q 2 ). (2) 

Therefore, in our notations 

gi(x,Q 2 ) = g q (x,Q 2 )+g g (x,Q 2 ). (3) 

There is no difference between our approach and DGLAP in this respect. The main difference between the IREE 
and the DGLAP equations is that due to the fc_L-ordering which is the key point of DGLAP, the ladder partons with 
minimal transverse momenta should always be in the lowest rung. However, such an assumption is an approximation 
which is true only for large x. When x is small, such partons can be in any rung in the ladder. More generally, 
each rung consists of two ladder partons which can be either quarks or gluons. Each of these two cases yields DL 
contributions and therefore have to be accounted for. It is convenient to write down the IREE in the to -space related 
to the momentum space through the asymptotic form of the Sommerfeld- Watson (SW) transform for the scattering 
amplitudes. However, the SW transform is defined for the signature amplitudes Mw. In particular for the negative 
signature amplitudes it reads (we drop the superscript "(-)" &s we do not discuss the positive signature amplitudes 
in the present paper): 



J., 

— {s/n 2 YZ{u)F r {u,Q 2 ) (4) 

-too 2m 

where r = q,g and £(w) is the negative signature factor, £(u>) = [1 — e~ I7rw ]/2 ss ittuj/2. It is necessary to note that 
the transform inverse to Eq. Q involves the imaginary parts of M r : 



2 f°° 

F r (uj,Q 2 ) = — / dpe~ puJ< 3M r (s,Q 2 ) 



(5) 



where we have used the logarithmic variable p = ln(s/^ 2 ). We will also use another logarithmic variable y = ln(Q 2 //x 2 ). 
In these terms, the system of IREE for F r (w, Q 2 ) can be written down as follows: 

(w + ^)F q (u, y) = -^2 [F qq {u)F q {u, y) + F qg (u))F g (uj, y)} , 

(w + ^)F g (u,v) = ^[F gq {u)F q {u,y)+F gg {u)F g {u,y)] (6) 
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where the anomalous dimensions Fik, with i, k = q, g, correspond to the forward amplitudes for quark and/or gluon 
scattering, having used the standard DGLAP notations for the subscripts. It is convenient to absorb factors 1/87T 2 in 
Eqs. © into the definition of new amplitudes Hik related to by 



H ik {u) = (l/87r 2 )^ fc H 



(7) 



The lhs of Eqs. © corresponds to apply the operator —{j, 2 (d/d[i 2 ) to M^. The first (second) term in the rhs of 
each equation (jBJ) corresponds to the case when the ladder rung with minimal transverse momentum is made of 
a quark (gluon) pair. The Born contribution does not appear in the rhs of Eqs. © because the Born amplitude 
j^Bom _ anc j j^Bom _ e 2 S/ /( s — Q 2 + ie) and therefore it vanishes when differentiating with respect to /x. The main 
difference between the system of the Altarelli-Parisi (AP) eqs. for the singlet g\ and Eqs. (JSJ) is that in the AP eqs. 
the amplitudes Hik contain NLO terms with two-loop accuracy, whereas Hik in Eqs. © account for all NLO terms 
in DLA. Then, these NLO terms cannot be included within the DGLAP approach whereas in our approach we take 



them into account in Sect, 
expressions for g q and g g : 



3, using /x 2 -evolution. Solving Eqs. 10 and using Eqs. (|1I2|I . we arrive at the following 



g q (x,Q 2 ) = 



2vr? 



{l/x) u C + {w)e u+v + C-{uj)e n - y 



(8) 



(x,Q 2 



— {1/xY C+{lo)- 
.™ Zm L 



2H n 



Factors C±(lu) will be specified in Sect. IV. We have denoted here 



2H„ 



(9) 



X — Hgg H qq , 



(10) 



R = (H a 



H qq ) 2 + 4:H qg H g 



and 



n+ = - 



H 



'i 'i 



H gg ± 



(Hqq-H gg ) 2 +AH nn H, 



(11) 



(12) 



Obviously, il + > 0_. Eq. JSJ includes the coefficient functions C±(u>) that should be specified. This can be done 
by different ways. Below we obtain C± in terms of the phenomenological quark and gluon inputs and then calculate 
the perturbative contributions to these inputs. However before doing so, we first calculate the anomalous dimensions 
H ik , (i,k = q,g). 

III. CALCULATING H ih 

As explicitly shown in eq. (|12fl . fl± are expressed in terms of the amplitudes Hik- As we are going to calculate 
f2± to all orders, we consequently have to know H^ to all orders in a s . To this aim we write and solve some IREE 
for Hik- The lhs of the new IREE again corresponds to differentiating the amplitudes with respect to — /x 2 (9/9/i 2 ) 
(cf Eqs. El . The rhs include, besides the ladder contributions involving Hik, their Born terms 



TjBorn 
H ik 



a-ik/u, 

and non-ladder contributions Vik, which we specify later. The system of eqs. reads 



(13) 



l0H qq = 

U)Hgg = 



L0H q 

ujH„ 



iqq 



V, 



gq 



V, 



gg 



H, 



H nn H n 



qq 1 ^^qg ±± gq^ 
" Hgg + HgqHqgi 

Hqg{H qq + Hgg), 
H gq (H qq + Hgg) . 



(14) 
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The non-ladder DL terms Vik appear in Eqs. I|14|l when the parton with minimal transverse momentum is a non- 
ladder gluon 1 . Such a gluon can be factorized, i.e. its propagator is attached to external lines. As the factorized 
gluon bears a colour quantum number, the remaining DL contribution -defined H c for the sake of simplicity - gets 
also a coloured content. When the factorized gluon propagates in the t -channel, H c belongs to the octet t-channel 
representation of SU(3), whereas all Hik belong to the singlet representation. Therefore in order to solve Eqs. <|14[) 
one has to calculate the octet amplitudes first. Fortunately, they can be approximated quite well by their Born values 
as they fall with energy very rapidly. With this approximation one obtains 



V lk = (15) 

where 

C N 
m qq = m gg = ~2N 2 , m qg = nf—, m gq = -NC F . (16) 

where rif is the number of the flavours. We assume n/ = 4. Furthermore D(lu) in Eq. 1)15(1 accounts for the running 
QCD effects for Vik- According to Ref. 0] it is given by 



1 r°° 

D H = 262 / dpe-^^{{p + rf)/r{) 

where 77 = h\(p 2 / K 2 QCD ) and b = (33 - 2n/)/127r. 

All cijfc in Eqs. I|14|l are proportional to a s . However, in Ref. it is suggested that in contrast to the DGLAP 
prescription a s = a s (k 2 _) (k± is the transverse momenta of the ladder partons), which holds universally in every rung 
of the ladder, the arguments of a s are different for different kinds of the rungs. In particular, the expressions for the 
quark-quark and the gluon-gluon rungs include a s depending on the time-like virtuality of the intermediate gluon, 
leading to the following expressions for a qq and a gg : 



. ( n -A- n) 2 A- tt 2 n A- n 



(17) 



2tt 
where 



A(w)C F 2A(u)N 

a gg = , (,!») 



A(w) 



dpe~ 



(P + v) 2 



(19) 



The 7r 2 -terms appear in Eq. H19[) because a qq and a gg involve a s with the time-like argument. On the other hand, 
in the expressions for the rungs mixing quarks and gluons, a s depends the space-like virtuality of the ladder gluons, 
which is approximated by an expression similar to (|19f) without the 7r 2 -terms. So, we arrive at 



n f A'{uj) _ A'(u)C f 



2tt 



with 



1 



dpe~ 



■ujp 



-V Jo (P + V) 2 - 

Once aik and Vik are determined, it is easy to obtain explicit expressions for Hik- Indeed, defining bik as 



(20) 



(21) 



bik — o-ik + Vik 



(22) 



1 We use the Feynman gauge through the paper. 
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and 



Y =-]J - 2 (V + M + V [(w2 ~ 2(6 ^ + h ^ ))2 ~ 4(6 ^ ~ h ^ ~ 16b M ) /2 - (23) 

we obtain 



= i(- + y + ^_^£),^ = i(c, + y-^_>) j (24 ) 

TT _ _^m_ tt _ bqg 

n 9i — y ' 99 — Y ' 

As a matter of fact, the expression ((23(1 for y appears in solving Eqs. 1(14(1 as one of the four solutions of a biquadratic 
algebraic equation. The signs before the roots in Eq. 1(23(1 are chosen so to get a matching between Hik = dih/u) + 
O(l/oj 2 ) when oj > 1 and the Born contributions of Eqs. I(2l)(l . Now all ingredients for f2± are specified and fl± can 
be obtained with numerical calculations. 

IV. INTERCEPT OF gi 

Eqs. (JSJ) and © give general expressions for g\. These expressions account for both DL contributions and running 
a s effects. The integrands in Eqs. I|8I9|) contain the coefficient functions and exponentials. In the present Sect, we 
consider the exponentials. We will study the coefficient functions in the next Sect. When the limit x — > is considered, 
one can neglect the exponents with f2_ in these equations and simplify the expression for Q + . Indeed, the behaviour 
of gi in this limit is driven by the leading singularity in Eq. 112|) . The singularities of fl + are related to the branching 
points of the square root. Using Eq. 1)23(1 one can see that the leading singularity of fl+(oj) is given by the rightmost 
root, uiq of the equation below: 



uj 4 - 4(6 m + b gg )u? + 16(b qq b gg - b qg b gq ) = . (25) 



Therefore, Eqs. I|8I9|) predict that 



91 



cM(i/zr°(Q 2 /M 2 r o/2 (26) 



when x — > 0, where loq is discussed below and the factor C(a->o) in Sect. V. Eq. 1(2 5(1 can be solved numerically. In 
our approach, ujq depends on r/ = ln(/i 2 / Aq C d ) . The result of the numerical calculation for u>o = u)o(rf) is represented 
by the curve 1 in Fig. 1. This curve first grows with 77, achieves a maximum where approximately luq = 0.86 and 
smoothly decreases for large 77. In other words we have obtained that the intercept depends strongly on the infrared 
cut-off /x for small values of /1 and smoothly thereafter. Quite a similar situation was occurring in Refs. [16J for the 
intercept of the non-singlet structure function g± ■ We suggest a possible explanation for this effect. The cut-off /x 
is defined as the starting point in the description of the perturbative evolution. Everything that affects the intercept 
at scales smaller than \i is attributed to non-leading effects and/or non-perturbative contributions. Had they been 
accounted for, the intercept would have been /i- independent. Without those non-leading/non-perturbativc effects 
taken into account, we then observe an important /i-dependence, which becomes weaker for large /i. The maximal 
perturbative contribution to the intercept A5 of 171 is obtained from the maximal value of ujq. Therefore we estimate 
the intercept as 

A s = max(wo(r})) ~ 0.86 . (27) 

Eq. ((27(1 includes the contributions of both virtual quarks and gluons. These contributions have opposite signs and 
partly cancel each other. It is interesting to note that when only virtual gluon contributions are taken into account, 
this purely gluonic intercept, A g is given by the maximum of the curve 2 in Fig. 1, which is slightly greater than 1. 
This value obviously exceeds the unitarity limit, similarly to the intercept of the LO BFKL|13|, though in a much 
softer way. Fortunately, by including also the contributions of the virtual quarks the intercept decreases down to A5 
of Eq. 1(27(1 . without violating unitarity. 
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We would like also to add the following observation. The total resummation of the DL contributions to g± performed 
in Refs. 0, was used in Refs. where a s was suggested to be fixed at the scale Q 2 , as in the framework 

of DGLAP. In other words, the parametrization a® L = a s (Q 2 ) was used in the DL expression for the intercept, 
Aj? L = 3.5y/3aP L /(2n) obtained in Ref. 0. Now, comparing this expression and Eq. 1|27[) . one easily concludes that 
the effective scale of a® L has no relation with Q 2 . This is in a complete accordance with the results of Refs. fl5|.[l6||. 
Indeed, the parametrisation a s — a s (Q 2 ) is valid only when the factorisation of the longitudinal and transverse 
momentum space is assumed. This approximation is correct for the kinematic Al region of large x and does not hold 
for small x. 



V. COEFFICIENT FUNCTIONS C*± 

In order to find g q and g g in Eqs. (JHHSJ), one has to finally specify the coefficient functions C±. One way to do it is 
to impose the matching 



9 q (x, Q 2 )Iq 2 = m 2 = A l( x o), 9g(x, Q 2 )|q2 =a1 2 = Ag(x ) 



(28) 



where xq — M 2 / s ! Aq(xg) and Ag(xo) are the initial densities of the polarized quarks and gluons respectively. Using 
the integral transform of Eqs. I|8I9[1 at Q 2 — /j 2 in the w-space the matching of Eq. (|28|l can be rewritten as 



C+ + C_ = Aq, a 



X + ^fR „ X-VE 



2H n 



2H„ 



= Ag, 



(29) 



L qg " L± qg 

with both Aq and Ag depending on u>. Combining Eqs. (JSJUJ) and (|29|l . we express C± through Aq, Ag and arrive at 



f lo ° rii i 



Ag + BAg]e 



A^Aq-BAg)e n - y ] , 



(30) 



2m 



EAq + A^Ag)e n+v 



- EAq + A^Ag 



(31) 



with 



±JL- 

2 2y/R 



g _ Hqg _ Hgq 



(32) 



Eqs. (|371)) and (fTTTf) describe the Q 2 evolution of the polarized quark and gluon densities g q and g g from Q 2 = [i 2 , 
where they are Aq(fj, 2 /s) and A<?(/i 2 /s) respectively, to the region Q 2 3> /i 2 . In the framework of our approach, /j, 2 
is the starting point for the Q 2 -evolution. Then one option would be to fix the initial parton densities Aq and Ag 
from phenomenological considerations. 

In this case, one can use Eqs. <|30I31I32[) in order to fix the sign of g\ in the small- x region. The problem of the sign 
involves an interplay of the gluon and quark contributions to <?i (see e.g. Refs. [^.pjf'l. Let us estimate the sign of g\ 
in the small- a; region, calculating the asymptotics of B, E. When x — ► 0, the main contributions in Eqs. (|30I31(I 

come from the terms proportional to exp[f2 + y]. According to the results of Sect. 4, the small- x asymptotics of g\ is 



g x ~ [{A { s ] + E s )Aq + (4 +) + B s )Ag] (l/x) A ° (Q 2 / > 2 ) A */ 2 



(33) 



where we have provided the factors A( ' , B, E with the subscript U S" in order to show their explicit dependence on 



lo = As- Substituting the numerical values A s 
arrive at 



(-) 



-0.31, A { + ] 



1.31, B s = 0.52, E s = -0.79 into Eq. we 



gi ~[-lAAq+1.8Ag] (l/x) As {Q 2 /^ s/2 



(34) 
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As Ag is positive, g\ is negative when 

Ag > 1.7 Ag . (35) 

Eqs. (|34I35I) are expressed in terms of the quark and gluon initial densities Ag, Ag defined (see Eq. I|28|)) at the 
scale x w /i 2 /s. On the contrary, the standard DGLAP expressions for gi involve the initial densities Sq and Sg 
defined at the scale x ~ 1. In order to compare our results to DGLAP, we should express Aq and Ag in terms of Sq 
and <5g. In the w-space, it means expressing Aq(u) and Ag(u>) through 8q{u>) and 8g(u>). We can do it within our 
framework by using the evolution of Ag and Ag with respect to s at fixed Q 2 w /i 2 from the starting point s w /z 2 to 
the region s>/i 2 . The system of the IREE for Ag and Ag is similar to Eq. © with the exeption of the fact that Ag 
and Ag do not depend on Q 2 and therefore the IREE for them do not include the derivatives d/dy. Therefore, we 
arrive at the following algebraic IREE: 

AgH = (e 2 /2)8q(io) + (1/lu) [H qq (uj)Aq(uj) + H qg {u)Ag{w)) , 

Ag{u) = (e 2 q /2)Sg(cu) + (l/tu)[H gq (uj)Aq(uj) + H gg (uj)Ag(tu)} . (36) 

where e 2 is the sum of the quark electric charges (e 2 = 10/9 for rif = 4), 8q is the sum of the initial quark and 

antiquark densities and Sg = — (A'{ui)/2-KLd 2 )Sg (A'(u>) is defined in Eq. JUJ) is the starting point of the evolution 
of the gluon density Sg. Eqs. (|36[) describe the s-evolution of the quark and gluon inputs from s ~ /i 2 to s » /j 2 . 
Solving Eqs. JH, 

we obtain: 

a [uj(uj-H gg )5q + ujH qg 8g] 

Ag = (e q 2)Y- — --, , (37) 

[w 2 - u{H qq + H gg ) + (H qq H gg - H qg H gq )\ 

9/ p _ W (ff g9 + + {Hqq H gg - H qg H gq )] ■ 

Substituting Eqs. (|37I38(I into Eq. I|34(l allows to express asymptotics of gi in terms of Sq, Sg: 

5 i - (l/2)[-1.2<5g - 0.08(5g ] (l/x) As (Q 2 / ^ 2 ) As/2 . (39) 

In contrast to Eq. (|34|l which involves the parton distributions Ag and Ag defined at the low- a: scale, at x ~ /i 2 /s, 
Eq. H39J) is expressed in terms of the densities <5g and Sg at x ~ 1. These densities can be fixed from phenomenological 
considerations. Eq. (|39|l states that gi can be positive at x <C 1 only when 5g is negative and large: 

\Uq +Sg < . (40) 



VI. SUMMARY AND OUTLOOK 



By constructing and solving the system © of infrared evolution equations, we have obtained the explicit expressions 
(|30I31|I for the polarized quark and gluon distributions g q (x,Q 2 ) and g g (x,Q 2 ). These expressions account for DL 
contributions to all orders and, at the same time, for the running a s effects. Both g q (x,Q 2 ) and g g (x,Q 2 ) depend 
on the low- a; quark and gluon distributions, Ag and Ag. Eq. (|34|) shows that g\ has the negative sign in the small- 
x region if Ag > X.lAg. By the same method, and using the s-evolution at Q 2 fts ^ 2 , the distributions Ag and 
Ag are expressed in Eq. I|37I38|I in terms of the initial densities Sq and Sg which are supposed to be fixed from 
phenomenological considerations at x ~ 1. Eq. I|39() shows that g\ is positive when the initial gluon density is negative 
and large: Sg < —158g, otherwise g\ is negative. 

We obtain that the expressions H30I31(I lead to the Regge behaviour (|34I39J) of g\ when x — > 0. The value of the 
intercept As depends on the infrared cut-off fi, but this dependence is quite weak for /i much larger than Aqcd- The 
estimated value of the intercept is given by Eq. I|27(l . and is a factor of 2.2 larger than the non-singlet intercept. The 
value of As = 0.86 is in a good agreement with the estimate As — 0.88 obtained in Ref. [2(j from analysis of the 
HERMES data. Our results also show that accounting for the gluon contributions only one would obtain a value of 
the intercept exceeding unity and therefore violating unitarity - similarly to the LO BFKL intercept - whereas the 
inclusion of the quark contributions stabilises the result. We prove that it is unrealistic to combine the resummation 
of the DL contributions to g\ with the DGLAP-like parametrisation for a s in the expressions for the intercepts. 

Finally, we note that it would be very interesting to implement our results by non-perturbative (lattice) calculations 
in order to check explicitly the independence of the total intercept on /i. 
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VIII. FIGURE CAPTIONS 

Fig. 1: Dependence on r\ of the rightmost root of Eq. I|25|l . Uq. Curve 2 corresponds to the case when gluon 
contributions only are taken into account; curve 1 is the result of accounting for both gluon and quark contributions. 
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